GRB 990712: First Indication of Polarization Variability in a 
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ABSTRACT 

We report the detection of significant polarization in the optical afterglow of GRB 990712 on 
three instances 0.44, 0.70 and 1.45 days after the gamma-ray burst, with (P, 6) being (2.9%±0.4%, 
121.1° ± 3.5°), (1.2% ± 0.4%, 116.2° ± 10.1°) and (2.2% ± 0.7%, 139.2° ± 10.4°) respectively. 
The polarization is intrinsic to the afterglow. The degree of polarization is not constant, and 
smallest at the second measurement. The polarization angle does not vary significantly during 
these observations. We find that none of the existing models predict such polarization variations 
at constant polarization angle, and suggest ways in which these models might be modified to 
accommodate the observed behavior of this afterglow. 
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The radiation from gamma-ray burst after- 
glows has been hypothesized from early on to be 
synchrotron emission from rclativistic electrons 
(e.g., Rees and Meszaros 1992; Meszaros and Rees 
1997). The observation of broad-band afterglow 
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spectra (Galama et al. 1998; Bloom et al. 1998) 
and the fact that the temporal evolution of the 
flux from radio to X-ray wavelengths often follows 
simple model predictions support this theoretical 
framework (Sari 2000; Piran 1999; van Paradijs 
et al. 2000). Since synchrotron radiation is in- 
trinsically highly polarized (with degrees of linear 
polarization up to 60-70%, Rybicki and Lightman 
1979), it would be natural to expect polarization 
in afterglow emission. 

The first attempt to measure optical polariza- 
tion of afterglow emission, for GRB 990123, set an 
upper limit of 2.3% (Hjorth et al. 1999). The next 
attempts to measure polarization were successful, 
leading to the detection of polarization of 1.7% in 
the afterglow of GRB 990510 (Covino et al. 1999; 
Wijers et al. 1999). Such low percentages of polar- 
ization are somewhat surprising, and models have 
been mostly concerned with trying to explain the 
depolarization of gamma-ray bursts. Two princi- 
pal hypotheses have been advanced. First, that 
the emission from the GRB is so highly ordered 
that the symmetry of the source causes the net 
polarization to average out to zero, even though 
it is locally high everywhere (Medvedev and Loeb 



1999). Second, that the emitting medium is very 
random on smaU scales, and that the emission we 
see is composed of many uncorrclatcd polarization 
patches so that the mean is again close to zero 
(Gruzinov and Waxman 1999). Both these models 
easily obtain the low polarization levels observed, 
and they differ only on the circumstances under 
which polarization may be observed, and what its 
temporal evolution should be. In the afterglow of 
GRB 990510, the close spacing of the early obser- 
vations and the large measurements errors later on 
precluded any strong conclusions about preferred 
models. 

GRB 990712 was detected with the BeppoSAX 
Gamma-Ray Burst Monitor and Wide Field Cam- 
eras Unit 2 on 1999 July 12.69655 UT (Heise et al. 
1999). The first optical follow-up observations 
started 4.16 hours after the burst, leading to the 
discovery of an optical transient (OT) of magni- 
tude Rot ~ 19.4 (Bakos et al. 1999). The light 
curve of the OT is quickly dominated by the light 
of its fairly bright host galaxy, with Rh ~ 22. 

A detailed description of the photometry of the 
OT is given by Sahu et al. (2000) and Hjorth et al. 
(2000). Galama et al. (1999) determined a redshift 
of z = 0.430 ± 0.005 for the OT. Spectroscopic 
observations of GRB 990712 are reported in Hjorth 
et al. (2000) and Vreeswijk et al. (2000). 

The outhne of this paper is as follows: in Sect. 2 
we describe the observations and data reduction. 
We analyze the results in Sect. 3 and discuss their 
significance for existing models in Sect. 4, as well 
as possible modifications to those models. Finally, 
we summarize our findings in Sect. 5. 

2. Observations 

Three epochs of observations were taken on 
July 13 and 14 with the Very Large Telescope 
(VLT) of the European Southern Observatory 
at Paranal, Chile, using VLT Unit Telescope 1 
{Antu) and the FOcal Reducer low dispersion 
Spectrograph (FORSl). All images were taken 
with a Bessel R band filter. In order to obtain the 
degree of linear polarization, a WoUaston prism 
and a half wavelength phase retarder plate were 
used as polarization optics. The WoUaston prism 
separates the incident light into two components 
(an ordinary and an extra-ordinary component), 
where the half wavelength plate is used to deter- 



mine which Stokes parameter is measured (U or 
Q). A mask producing 22" wide parallel strips 
was used to avoid overlap of the ordinary and 
extra-ordinary rays. Each observation consisted 
of four exposures centered on the position of the 
OT, with the phase retarder plate rotated by 22.5° 
between successive exposures. 

The data were reduced in a standard way 
with the NOAO IRAF^ package ccdred, first bias- 
subtracted and then fiat-fielded. The fiuxes of the 
point sources in the field were determined using 
aperture photometry, with an aperture radius of 
one FWHM of the Point Spread Function. 

To correct for any instrumental or local inter- 
stellar polarization, we measure the polarization 
of the OT relative to 21 field stars in the same 
image. We have plotted the Stokes parameters U 
and Q of both the field stars and the OT in Figure 
1. The OT clearly stands out, having the lowest 
U value. We verified that there are no systematic 
variations of U and Q of the field stars with mag- 
nitude or position on the CCD, and therefore we 
can correct the polarization parameters of the OT 
for foreground effects by subtracting the mean U 
and Q values of the field stars from those of the 
OT. 

3. Results 

To derive the degree of linear polarization, we 
first calculate Q and U using standard equations 
(see, e.g., Ramaprakash 1998), both for the OT 
and for the field stars. We then subtract the av- 
erage Q and U values of the field stars from those 
of the OT. The resulting Q and U are used to 
calculate the degree and the position angle of the 
linear polarization. We have also used a different 
method (see, e.g., di Serego Afighieri 1997), which 
within the errors leads to the same results. This 
method uses the relation S{<j)) = P cos 2(0 — (j)), 
with the parameter S{(f>) being a measure for the 
ratio between the two components of the incident 
light, separated by the WoUaston prism, and the 
corresponding angle of the prism. Figure 2 gives 
a cosine fit to the data of the first epoch. 

In Figure 3 we have plotted the polarized flux 
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together with the R- and V-band hght curve. The 
plot clearly shows the change in polarization, while 
the light curve exhibits a smooth decline. There is 
also no indication of a change in the (V-R) color 
of the afterglow. 

We assume that the degree of linear polariza- 
tion and the position angle are constant during the 
observation. The total intensity of the transient is 
not, as it is decaying according to a power-law, 
/ ex t-", with a = 0.97 (Sahu et al. 2000). The 
fact that the powerlaw index is known gives us 
the possibility to correct for the decaying inten- 
sity. These corrections turn out to be so small 
that we don't find any differences in the polariza- 
tion between an assumed constant intensity from 
the start to the end of the observation, and a 
powerlaw-like declining intensity. 

Both the degree of linear polarization and the 
position angle versus the time since the burst are 
plotted in Figure 4. We see that the polariza- 
tion percentage, P, decreases between 0.44 and 0.7 
days after the burst, with 3.2cr significance. Then 
at 1.45 days, P is greater again, but since the dif- 
ference between the last and middle observation is 
only 1.5(7, we cannot be too sure about this rise of 
P. 

In case that the degree of polarization would 
be constant, with a mean value of 2.1% ± 0.3%, 
the P value at epoch 1 is just within Str above the 
mean, where for epoch 2 it is just 3cr below it. The 
variability could then be caused by a systematic 
error. To check this, we have compared the polar- 
ization values of the separate field stars at each of 
the three epochs with the mean value. We have 
not found any evidence for such a systematic error, 
and conclude that the observed variability in the 
polarization is most likely intrinsic to the source. 

Therefore, we think that we have clearly de- 
tected, for the first time, variation in the polariza- 
tion of a GRB afterglow, which could either be a 
decline that is initially steep and then levels off, 
or a decline followed by a rise. The polarization 
angle, at the same time, never changed by more 
than one standard deviation from one observation 
to the next. 

As argued for the case of GRB 990510, polar- 
ization in a rapidly varying source is unlikely to 
be induced by interstellar scattering (Wijers et al. 
1999). In the present case, where the polarization 



itself varies, this is even more strongly so: stars 
with polarization induced by interstellar scatter- 
ing are favorite polarization standards, because 
their degree of polarization is very constant. 

4. Discussion 

Intrinsic polarization from a synchrotron source 
can be as large as Pmax ~ 60-70% (Rybicki and 
Lightman 1979) if the magnetic field is oriented in 
one direction. However, for an unresolved source 
the net polarization will be small if the different 
directions of the polarization average out. This 
could be caused by highly tangled magnetic fields, 
or by a very simple symmetry in the large-scale 
field pattern. 

If the magnetic field in the GRB afterglow is 
highly tangled, we can think of it as a source 
consisting of N patches (Gruzinov and Waxman 
1999). The net polarization resulting from the 
source will then be of order Pmax/"sfN ■ The max- 
imum degree of linear polarization observed (at 
epoch one) requires a magnetic field divided up 
into ~ 400 patches. For the second epoch, we see a 
decrease in the degree of linear polarization (with 
9 being constant), which would lead to ^ 2500 
patches, six hours later. We note, however, that 
in that case the polarization angle should also vary 
by the same percentage, implying that we expect 
a change of order a radian between the first and 
second epoch in a tangled-field model. However, 
we see no significant change in the polarization an- 
gle, with a 1-radian variation ruled out at the Aa 
level, making this model unlikely. 

A symmetric model for the polarization can- 
cellation arises when the magnetic field is ordered 
perpendicular to or parallel to the ring of emission 
we see from the afterglow at any given time (Wax- 
man 1997; Panaitescu and Meszaros 1998; Sari 
1998). This could naturally arise in some insta- 
bilities that generate magnetic fields (Medvedev 
and Loeb 1999; Gruzinov 1999). In a spherically 
symmetric fireball, the polarization would then be 
exactly zero, so an extra effect is needed to break 
the symmetry and get a net polarization. One 
possibility is that some turbulence induces bright- 
ness variations, thus weighting some polarization 
directions more, or that an external effect such 
as scintillation or microlensing might enhance the 
emission from some parts of the ring (Loeb and 



Perna 1998; Medvedev and Loeb 1999). However, 
any polarization variations from such a mechanism 
would be expected to be both in degree and angle, 
contrary to what we see. 

Another symmetric model is a jet: as explained 
in Ghisellini and Lazzati (1999) and Sari (1999), 
a coUimated burst will naturally exhibit polariza- 
tion, up to 20%, around the time when the light 
curve steepens. However, we have no evidence of 
such a break in the light curve of GRB 990712. 
Also, in a jet model, the degree of polarization 
varies without change of the polarization angle, 
until the polarization goes through zero and the 
angle suddenly changes by 90 degrees. Therefore, 
one could in principle have a situation of vary- 
ing polarization without variation of the angle as 
we see here in a jet model. However, the sharp 
drop from epoch 1 to 2 followed by little change 
to epoch 3 seems hard to get without a zero tran- 
sit around epoch 2 if we compare Figure 4 with 
the theoretical curves of Sari (1999, Fig. 4) and 
Ghisellini and Lazzati (1999, Fig. 4). So at least 
for the simple models published thus far, our data 
cannot be explained by a beamed jet. 

In summary, current models of polarization 
variations predict significant changes of the polar- 
ization angle along with variations of the degree 
of polarization, or the degree of polarization to be 
around a peak for a constant polarization angle in 
case of a jet model. Since we do not observe such a 
combination of the polarization angle and degree 
of polarization in the afterglow of GRB 990712, 
we conclude that no current model adequately ex- 
plains our data. 

5. Conclusions 

We have observed significant polarization for 
the afterglow of GRB 990712 during three epochs 
from 0.4 to 1.5 days after the burst. The polar- 
ization percentage varies by 3.2a from 3.0% at 
the first epoch to 1.2% at the second, while the 
polarization angle remains essentially constant. 
We show that neither tangled-field nor broken- 
symmetry models of the polarization of afterglows 
can explain this constancy of the polarization an- 
gle while the degree of polarization goes up and 
down. It appears that the afterglow polarization 
has some amount of memory for direction, while 
varying in strength. One might speculate that 



this is telling us something about the formation 
of the magnetic field in the shocked afterglow ma- 
terial. Possibly, this field is grown from a seed 
field that is embedded in the swept-up ambient 
material. While the amount of amplification and 
field strength may vary with time, and give rise 
to variable polarization, the direction of the seed 
may remain imprinted on the amplified field, and 
thus the polarization angle may be relatively con- 
stant. (Few-day old afterglows are 10^^ — 10^^ cm 
in size, on which scale interstellar magnetic fields 
can be coherent.) 

The results show that in order to further test 
models, we need more measurements per burst, to 
better sample the polarization behavior, and pos- 
sibly over a larger range of time. Since the mea- 
surements require photometry with a S/N of 300 
or so, they already stretch the capabilities of the 
VLT after 1.5-2 days, so it is unlikely that we shall 
be able to do measurements of the polarization at 
later times. However, it may be possible to extend 
the time interval of polarization measurements to 
earlier times. In a power-law process like GRB af- 
terglows, it may be as profitable to move the first 
observing time forward from 0.4 to 0.2 days as it 
is to extend the last one from 2 to 4 days. 
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Fig. 1. — (Q, U) plot of the stars in the field after 
instrumental correction, for epoch 1, 2, and 3 (top, 
middle and bottom respectively). The OT clearly 
stands out with respect to the other stars in the 
field, indicating that it has a considerable intrin- 
sic amount of polarization. {Q, U) averages of the 
field stars for the three epochs are (0.0194±0.0004, 
-0.0021 ± 0.0005), (0.0200 ± 0.0003, -0.0024 ± 
0.0004), (0.0190 ± 0.0005, -0.0019 ± 0.0004) re- 
spectively, showing that there is no significant 
change in the interstellar and instrumental polar- 
ization. 



0.04 



0.02 






■0.02 



-0.04 



"I 1 1 1 1 1 1 1 \ r 

P = ( 2.3 ± 0.7 )% 
0= ( 148.2 ± 9.5 Y 




J I I L 



J I \ L 



J I I L 







50 100 

0(degrees' 



150 



Fig. 2. — The parameter S{(f) at four different 
position angles (p. The data are fit with a cosine 
function. The amphtude of the fit corresponds to 
the degree of linear polarization, and its niaxinium 
gives the position angle of the polarization. 
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Fig. 3. — The light curve of GRB 990712. V and R band data (filled circles and triangles respectively) are 
from Sahu et al. (2000). The open circles indicate the polarization magnitudes. 



Table 1 
Log of the polarimetric observations of the GRB 990712 afterglow with ESO VLT Antu 
(UTl) AND FORSl. The exposure time for epoch 1 and 2 was 4 x 300 seconds and for the 

3rd epoch 4 X 450 seconds. Pav and Oav DENOTE THE AVERAGE VALUES OF THE FIELD STARS (nOT 
CORRECTED FOR INSTRUMENTAL AND INTERSTELLAR POLARIZATION). 



UT date^ 
(1999) 



Days 
after the burst 



seeing 
(arcscc) 



P 

(percentage) 



(deg) 



(percentage) 



(deg) 



July 13.1402 
July 13.3936 
July 14.1415 



0.4436 
0.6971 
1.4450 



1.0 
0.8 
1.1 



2.9 ±0.4 
1.2 ±0.4 
2.2 ±0.7 



121.1 ±3.5 
116.2 ±10.1 
139.2 ±10.4 



1.89 ±0.02 
2.06 ±0.02 
1.91 ±0.02 



175.5 ±0.3 

176.3 ±0.2 

175.4 ±0.2 



^Mid-exposure date. 

''The seeing varied somewhat during the measurements; the values listed here are averages. 
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Fig. 4. — The variation of P and 
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during the three observations. 
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